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Abstract 

Constructing the reduced density matrix for a system of three massive spin— \ particles de- 
scribed by a wave packet with Gaussian momentum distribution and a spin part in the form of 
GHZ or W state, the fidelity for the spin part of the system is investigated from the viewpoint 
of moving observers in the jargon of special relativity. Using a numerical approach, it turns out 
that by increasing the boost speed, the spin fidelity decreases and reaches to a non-zero asymptotic 
value that depends on the momentum distribution and the amount of momentum entanglement. 

keywords:Wigner rotation, spin density matrix, Gaussian momentum distribution, fidelity, GHZ 
state, W state. 



1 Introduction 

The role of special relativity in framing statements about quantum information is illustrated by the 
fact that quantum entanglement can depend on the reference frame of the observer. In practice, 
Lorentz transformations can change the entanglement of the spins of massive particles. Relativistic 
effects on quantum entanglement and quantum information is investigated by many authors. One 
of the early works in this area has considered a single free spin-^ particle and by calculating the 
reduced density matrix, it is shown that the spin entropy is not a relativistic scalar [1J. Alsing and 
Milborn [2J studied the Lorentz transformation of maximally entangled Bell states. They concluded 
that entanglement is Lorentz invariant. The entanglement between the spins of a pair of particles may 
change because the spin and momentum become mixed when viewed by a moving observer [3j. Li an Du 
have investigated the quantum entanglement between the spins of spin-^ massive particles in moving 
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frames, for the case that the momenta of the particles are entangled |4j. They have shown that, if the 
momenta of the pair are appropriately entangled, the entanglement between the spins of the Bell states 
remains maximal when viewed from any Lorentz-transformed frame. Bartlett and Terno showed that 
relativistically invariant quantum information can be encoded into states of indistinguishable particles 
[5] . Recently, simple examples have been presented of Lorentz transformation that entangle the spins 
and momenta of two spin-^ particles with positive mass such that no sum of entanglements have been 
found to be unchanged [BJ. Fidelity for the spin part of a system of two spin— | particles described 
by a Gaussian momentum distributed wave packet is studied from the view point of moving observers 
and it is shown that the fidelity decreases by increasing the boost velocity [7j. Bell's inequality in 
moving frames has been considered in several papers [HJ [91 [TOl [TTl [12l I13| . The degree of violation 
of Bell's inequality will decrease with increasing the velocity of the observers if the directions of the 
measurements are fixed. However, this doesn't imply a breakdown of nonlocal correlation since the 
perfect anti-correlation is maintained in the appropriately chosen different directions. Some efforts 
have been done for extending these ideas to tripartite systems. For example, Lorentz transformation 
of three-qubit GHZ state is studied and it is shown that Bell's inequality is maximally violated for 
this state [2]. In tripartite discrete systems, two classes of genuine tripartite entanglement have been 
discovered, namely, the Greenberger-Horne-Zeilinger (GHZ) class [151 05] an d the W class [TTl HH] - 
Some authors have provided proposals for generation and observation of GHZ or W type entanglements 
[H [20 [m [221 [23]. 

In this paper we consider a moving system containing three spin-| massive particles such that in 
the rest frame, its spin part be entangled as one of the GHZ or W states. In the present approach, 
we introduce a Gaussian momentum distributed wave packet represented in the momentum space 
for the system as viewed in the rest frame. Also we introduce the wave packet for the system as 
viewed by a boosted observer. Then we focus on the spin part of each wave packet by finding the 
corresponding reduced density matrix. As a result of relativistic spin decoherence, the reduced density 
matrix observed in the boosted frame is mixed even though it is prepared to be pure in the rest frame. 
We quantify the amount of mixing via calculating the fidelity for these two reduced density matrices. 
We will consider general boosts in the xz-plane. To see the effect of momentum entanglement on the 
results, the momentum part is chosen to be in extreme cases of momentum product or momentum 
perfectly correlated. 



2 Lorentz transformation of reduced density matrix 

A three-particle quantum state is expressed by 

|Pi,o-i;P2>°'2;P3 ) '3) = a t (p 1 ,cri)a t (p 2 ,CT2)a t (p3,o- 3 )|$o), (1) 

where p is the 3- momentum vector, a is the spin label, is creation operator and |$o) is the Lorentz 
invariant vacuum state. The state ([I]) has a Lorentz transformation property |24| as 

U(A)\ Pl ,a r ,p 2 ,a 2 ;p 3 ,a 3 ) = £ D a , ai (W(A, Pl ))D a , 2a2 (W(A, P2 ))D a/sa3 (W(A, P3 )) (2) 

x\Ap 1 ,a' l ;Ap 2 ,a' 2 ;Ap 3 ,a 3 ), 
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where Ap is the spatial part of Ap, D a / a (W(A,p)) is the unitary representation of the Wigner rotation 
operator, and W(A,p) is the Wigner's little group element 

W(A,p) = L- 1 (Ap)AL(p), (3) 

where L(p) is the standard boost that takes a massive particle of mass m from rest to a 4-momentum 
p. The transformation of the creation operator is as 

L/(A)at(p, f r)[/- 1 (A) = £ D a , a (W (A, p))a^ (Ap, a'). (4) 

a' 



Let p be the unit vector along the 3-momentum of a particle as viewed in the rest frame, and 
consider a boost along e with speed V. Then, the Wigner rotation operator is represented as 

D(W(A,p)) = Icos^ + i(a • n)sin^, (5) 

where 1 is the unit 2x2 matrix, a = (a\, 02, 03) denotes the Pauli matrices, n = ig x & and 

Q coth $ coth ^ + e • p 

cot - = ^-77 h -, 6) 

2 |e x p| w 

where 

p° V 
cosh£ = — , tanhr/ = — . (7) 

m c 

In the rest frame of the observer, the 4-momentum of each particle can be written in polar coor- 
dinates as 

p^ = [mc cosh £, mc sinh £(cos 93 sin 1?, sin (p sin cos 1?)] , (8) 

where tanh^ = - and v is the speed of the particle. To be specific, we suppose that particles are 
moving along the positive x-axis, i.e., ^ = § and ip = 0, then the 4-momentum reduces to 

p^ = (mc cosh £, mc sinh £, 0, 0). (9) 

This means that all of the three particles are assumed move along the positive x-axis, so for an arbitrary 
boost direction e = (cos 4> sin 9, sin <p sin 9, cos 9), the axis of Wigner rotation n is perpendicular to the 
x-axis in the same one direction for all three particles. We will consider only (p = 0. 

In this case the boost direction lies in the xz-plane and the Wigner rotation is about the y-axis as 
D(W(A, P)) = ( ftW %W ) = ( -| 2l ) < 10 > 



D it (Q) D u (n) J V -sm 



where 

17 coth $ coth ^ + sin 9 , 

cot- = ? 1 , (11) 

2 cosfl v ; 

One may provide the argument using 3-particle states ([1]) and ([2]) which require a sharp momentum 
distribution around momentum p for each particle. But the realistic situation involves a wave packet 
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of the system with a definite momentum distribution. We follow the argument in terms of Gaussian 
momentum distributed wave packets. 

In the rest frame of an observer the wave packet in momentum representation for the 3-particle 
system generally can be expressed as 



E 




d 3 p 1 d 3 p 2 (fp 3 g ai a 2 a 3 (Pl, P 2 , P 3 )|Pl; P2: °2! P 3 > °3>, 



(12) 



where <7o-io-2o- 3 (Pi) P 2l P3) is a distribution function for momentum and spin that is normalized as 



E 




(i 3 Pl (i 3 p 2 (i 3 P3|5 (Tl(T2(T 3(p 1 ,P2,P3)| 2 = 1, 



(13) 



which makes to have 



(Pi,o-' 1 ;p2,o-2;P3,cr3|p 1 ,(Ti;p 2 ,cj2;P3,cr3) = 5 3 (pi - Pi)5 3 (p' 2 - P 2 )£ 3 (p' 3 - P3)^A^ 2 <^ CT3 - (14) 
Now, regarding ([2]), for a boosted observer the state (fT2|) changes to 



1^ 



^ / / * /(AP1)° /(Aj*)° (Aps) 



EE/ d3 pi / d3 P2 / rf3 P31 



01CT203 CT ' (j'(jL 



n 



50-10-20-3 (Pl 5 P2 5 P3) ( 15 ) 



xD CT ^ 1 (W(A,p 1 ))D CT / CT2 (W(A,p 2 ))^ 3CT 3(W(A J p3))|Api^i;Ap 2 ,a 2 ;Ap3 J a 3 ), 
where D^^. is given by (fTUj) . 



The density operators corresponding to these pure states are p = and p b = \i(> b )(ip b \, however 

we need the reduced density operators g and g b obtained by tracing over the momentum of the density 
operators, that is 



d 3 p 1 d 3 p 2 d 3 p 3 g a > x a ^ 3 (p 1; p 2 , p 3 )g* ai a2(T3 (pi, p 2 , P 3 ), 



(16) 



and 

h 



Q^, m = E E / / / ^Pi^Pa^Ps (17) 

°1 CT 2°3 °1 2 3 

x C CT1CTi «(W(Ai,pi))D CT2CT »/(W(A^ 

Clearly g b denotes the Lorentz transformed form of g. It must be noted that g b will be mixed even if 
g is pure. 

Our purpose here is to calculate the fidelity of g and g b , which accordingly we call it the spin 
fidelity F s . Then, we use the Uhlmann definition |25[ 126] for fidelity of mixed states, that is 

(18) 



Tr(y -^fg g \fg) 



Fidelity is a basic ingredient in communication theory and for any given communication scheme it is a 
quantitative measure of the accuracy of the transmission. It takes numbers between and 1; a perfect 
communication corresponds to the fidelity 1. In the present argument the spin fidelity quantify how 
g b looks like g. 

Recall that in our problem each 3-momentum has only one component along the x-axis. Then all 
the three-dimensional integrals reduce to one-dimensional integrals where the bold notation for the 
momenta is suppressed. In the following p stands for the rr-component of 3-momentum, not for the 
4- momentum. 

Now we choose the distribution function g ai a2a 3 (pi,P2,P3) such that the spin and the momentum 
parts be separable and assume that the spin part is in GHZ or W state. In tripartite discrete systems, 
two classes of genuine tripartite entanglement have been discovered, namely, the GHZ class |15|, [T6] 
and the W class |171 118j . These two different types of entanglement are not equivalent and cannot 
be converted to each other by local unitary operations combined with classical communication. In 
terms of the spin basis, the GHZ state has the form |GHZ) = (| ft!) + I III)) an d the W state 
takes the form |W) = (| tH) + | 4-TI) + I lit))- The entanglement in the W state is robust against 
the loss of one qubit, while the GHZ state is reduced to a product of two qubits. According to the 
geometric measure of entanglement, the W state has higher entanglement than the GHZ state does 
[27] . Methods are proposed for generation and observation of GHZ or W type entanglements \22\ 123]. 
In our argument, it becomes apparent that there is an interesting contrast between the behavior of 
spin fidelity for these two states. 



3 GHZ state 



In this section we specify ga 1 a 2 o- 3 {pi,P2,P3) such that in the rest frame the spin part of the wave 
packet be in the GHZ state and, to perceive the effect of momentum correlation on the results, the 
momentum part be in states with different momentum correlation. 

First we choose the distribution function as 

9a 1 a 2 a i (Pl,P2,P3) = f (Pi ) / ) / (>3 ) (<W<W<W + <W<W<W) ' ( 19 ) 

which as applied in (PT21) . evidently gives a state that its spin part is entangled as the GHZ state and 
the momentum part is separable, i.e., the momentum entanglement is zero. We refer to this choice 
as the momentum product case. For evaluating the integrals, we should pick out a specified form for 
f(p). Here we consider it as 

(20) 

which shows a Gaussian distribution (minimum uncertainty) of momentum around p = with a width 



(a 2 7r) 



-exp 



1 (p 
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determined by a. Applying (fT9|) in (fT6|) and (fT7|) we get 



which is pure and 



Sa[ a' 2 a' 3 , o\ 02 <T3 



/10000001\ 
00000000 
00000000 

1 00000000 

2 00000000 
00000000 
00000000 

\l 1/ 



o / dpi|/(pi)| 2 / ^ 2 |/fe)| 2 / dp 3 |/(p 3 )| 



+^ it (fii)D ff , t (n 2 )£) o4t (n3)i>^(ni)^^n2)^(n3) 

+^(J2 1 )D CTU (0 2 )^(0 3 )iP* lt (J2i)C: 2t (0 2 ) J D* 3t (03) 
+^(J2 1 )L> ffU (fi 2 )L) CTU (03)^(J2i)L>: 2; (fi2)£';3;(03)], 

which is mixed. Since g is pure g 2 = g and (|18p can be written as 

2 



Tr(V gg b g) 



Using ([2T]) and ([22]) in ([23]) . the spin fidelity is obtained as 

1 



1 



We apply the Wigner rotation (jlOp in (|22p . then we obtain 



[ COS Qi COS fi 2 COS ^3 + cos fii cos Jl 2 + cos Ol COS 



where for each particle 

2 

cos fl = 



TV 71 " Jo 



die 



+COS fi 2 COS fls + COS fix + cos ^2 + COS + 1] , 

cos 2 6 (1 + cosh Tj cosh £ — cosh rj — cosh £) 



-7 sin. 



inh 2 £ 



coshcf 



1 



1 + cosh 7/ cosh £ + sin # sinh i] sinh £ 



(21) 



(22) 



(23) 



(24) 



(25) 



where 7 = a/mc, sinh£ = p/mc and we have used (fTTj) and (|20|) . Note that cosfii = cos$l 2 
then spin fidelity ([25]) reduces to 

1 



— ( cos ft + 3cos + 3cos $7 + 1 



(26) 



COS $^3, 



(27) 



Next, we consider a case that two of three momenta, say P2 and p 3 are perfectly correlated. We 
refer to it as the 2-momentum correlated case. Thus we write 

I 

9<na 2 <73(Pl,P2,P3) = f(Pl)f(P2)V S (P2 -P3)—7= (<W<W<W + <W<W<W) , (28) 
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which as substituted in (fT6|) and (fTT|) gives a pure g as ([2T|) and a mixed Applying the results in 
(|23|) . we have 



1 



-F s = — I COS $ll COS 2 $^2 + 2 COS ill COS f&2 + COS 2 0,2 + COS ill + 2 COS 0% + 1 1 , 



(29) 



which is comparable with (|25p . Again note that cos 11 1 = cosf^ then 



1 



F s = — ( cos 11 cos 2 H + 2 cos H + cos 2 + 3cos + 1 



(30) 



where 



cos 



2 



7V 71 " Jo 



^ e - 7 2 sinh 2 ^ cosh ^ 



cos 2 0(1 + cosh rj cosh £ — cosh r] — cosh £) 
1 + cosh 77 cosh £ + sin sinh rj sinh £ 



(31) 



Here we note that, in (|28p as well as in the following arguments, the delta functions should be 
regarded as limits of analytical functions under certain conditions. Precisely, perfectly correlated 
momenta should be regarded as a limiting case of entangled Gaussian momenta [3]. An experimental 
situation for generating the momentum entanglement is discussed by Lamata et.al. [28]. They studied 
the dynamics of momentum entanglement generated in the lowest-order QED interaction between two 
massive spin-^ charged particles, which grows in time as the two fermions exchange virtual photons. In 
this scheme the degree of generated entanglement between interacting particles with initial well-defined 
momentum can be infinite, however, they explained this divergence in the context of entanglement 
theory for continuous variables, and showed how to circumvent this apparent paradox. 

In order to have all the three momenta perfectly correlated ( 3-momentum correlated case) we 
choose 

0o-icr 2 <T 3 (Pl>P2,P3) = f(Pl)V S (Pl ~ P2)S(P1 ~ P3) (32) 

which as substituted in (fTBj) and (fT7]) . gives 

F s = - (cos 3 + 3cos 2 + 3cos~n + l) , (33) 

where 

o poo 

cos 3 U = — -= / di e-T" 2 sinh ? cosh i 
7V 7 *" Jo 



1 - 



cos 9 (1 + cosh 7] cosh £ — cosh rj — cosh £) 
1 + cosh r] cosh £ + sin 9 sinh rj sinh £ 



(34) 



There is no analytical solution for the integrals (|26p . (|3ip and (|34p hence we switch to a numerical 
approach for evaluating the spin fidelity F s . It reveals that the behavior of F s in terms of the boost 
velocity is in the same form for all the #'s and the most change occurs for 9 = 0, that is, when the 
boost is along the z-axis. Therefore, we proceed with substituting 9 = in the integrals. Fig. 1 
shows F s plotted numerically in terms of the boost parameter rj for a given width for the momentum 
distribution. The solid curve shows F s given by (f2"T|) in the momentum product case. The dashed 
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Figure 1: Spin fidelity F s versus the boost parameter rj in the GHZ case plotted for 7 = 20. The 
solid curve is plotted for the momentum product (zero momentum entanglement) while the dashed 
and the dashed-dotted correspond to the 2-momentum correlated and 3-momentum correlated cases, 
respectively. For small r\ the curves coincide however at large r] they slightly split. 



curve is plotted for F s given by (|29p in the 2-momentum correlated case and the dashed-dotted curve 
is for F s given by (133|) in the 3-momentum correlated case. We see that by increasing the boost 
velocity more spin decoherence occurs and expectedly the spin fidelity decreases with increasing r\. 
By increasing the momentum correlation, F s decreases less, such that for small rj the curves coincide, 
and for 77—7-00 (ultra relativistic limit) they slightly split to non-zero asymptotic values. It can be 
shown that by decreasing the width 7, the spin fidelity becomes less sensitive to 77, hence the slope of 
the curves decreases. 



4 W state 



As the previous section we investigate the W state in three cases of different momentum correlation. 
We begin with the momentum product case, by choosing 

9cT l a 2 <7 3 ipi,P2,Pz) = /(Pl)/(P2)/(>3)-/= (<W<W<W + <W<W<W + <W<W<W) • ( 35 ) 
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By this choice (fT6|) leads to 



6 = n 



/00000000\ 
00000000 
00000000 
1 1 1 
00000000 
1 1 1 
1 1 1 

Voooooooo/ 



(36) 



which is pure and (1171) gives 

tkoMw = \j *i|/(px)| a / d P 2\f(P2)\ 2 J d P3 \f( P3 )f (37) 

x[£> <t (ni)D^(n 2 )D^(n 3 ) + d^^d^^d^^) + d^^d^^d^^)} 

Again apply these operators in (f23|) and after doing some manipulation the spin fidelity is calculated 
as 



which as evaluated by ([TO]), becomes 



which reduces to 



Then, let 



— [21 COS ill cos cos r?3 + 5 COS ft\ COS ^2 + 5 COS ft\ COS $^3 

I A- 



+5 cos r^2 cos $^3 + 5 cos ill + 5 cos O2 + 5 cos SI 3 + 21] , 

1 / 3 2 \ 

— 7 cos Jl +5 cos £1 +5 cos f2 + 7 . 
24 V / 



Fs " 24 



(38) 
(39) 

(40) 



9a 1 a 2 a 3 (Pl,P2,P3) = f(p\)f(pl) \/S(j) 2 ~ Ps) (<W<W<W + <W<W<W + <W<W<W) ' ( 41 ) 

which apparently describes the 2-momentum correlated case. After doing some manipulations this 
leads to the following expression for the spin fidelity 



— (39 msQ cos 2 Q + 7 cos- O + [() Co s Q ~ - 3 cos Q + !<) 



(42) 



Finally, we consider the 3-momentum correlated case designated by 

Action 0-3 (Pl,P2,P3) 



P3) 



(43) 
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Figure 2: F s versus i] in the W case plotted for 7 = 20, for the boosts along the z-axis. The solid 
curve is plotted for the momentum product (zero momentum entanglement) while the dashed curve 
and the dashed-dotted curve correspond to the 2-momentum correlated and 3-momentum correlated 
cases, respectively. 



Using this, we find the spin fidelity as 

F s = (75 C os 3 ft + 25 cos 2 ft - 39 ooTft + ll) . (44) 



Again the most change in these fidelities occurs for = in the integrals and we present the 
result for this case. The curves in Fig. 2 are plotted numerically and describe the behavior of F s in 
the present case in terms of the boost parameter r\ for a given width 7. As is indicated, the solid 
curve, the dashed curve and the dashed-dotted curve are plotted for the momentum product case, 
the 2-momentum correlated case and the 3-momentum correlated case, respectively. Comparing with 
Fig. 1, we conclude that the F s again descends to nonzero asymptotic values but there is a significant 
separation between the three curves. This means that the present W case is more sensitive to the 
momentum entanglement. Also, note that the order of curves in Fig. 2 is inverse of the order of curves 
in Fig. 1. 
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5 Conclusions 



In this work we investigated a system of three massive particles described by a Gaussian momentum 
distributed wave packet such that in the rest frame its spin part was entangled as the GHZ state or 
the W state. Then we constructed the wave packet of the system as viewed by a boosted observer, 
by using the corresponding Wigner rotation operators. We focused on the spin part of the system by 
tracing out the momentum part and finding the reduced density operators both for the rest observer 
and the boosted observer. Using these reduced density matrices and the Uhlmann formula for fidelity, 
the spin fidelities were formulated separately when there was no momentum correlation, when two 
of three momenta were correlated and when all the three momenta were correlated. We could not 
evaluate fidelities analytically, so we utilized a numerical approach to plot F s in terms of the boost 
parameter rj, as Fig.s 1 and 2 show. 

We conclude that for the GHZ case, by increasing the boost velocity, F s falls to non-zero asymptotic 
values that increase as the momenta become more entangled. One may explain this behavior by 
regarding the results of the refs. [3] and [I]. By boosting the wave packet, we move some of the 
spin entanglement to the momentum part and simultaneously the momentum entanglement appears 
to to be moved to the spins. The amount of transferred entanglement grows with increasing the boost 
velocity. Tracing out the momentum from the Lorentz-transformed density matrix destroys some of 
the entanglement. This process causes to decrease the spin entanglement in the boosted frame. When 
the momenta are correlated, the transfer of momentum entanglement to spins compensates somewhat 
the decrease of spin entanglement and then the spin fidelity decreases less. However, as the figures 
show, this becomes more significant at large boost velocities. For the W case, the situation is inverse 
and by increasing the momentum entanglement, the spin fidelity decreases more. 
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